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ASYMPTOTIC  SOLUTIONS  FOR  THE  EXPANSION  OF  GAS  INTO  VACUUM 

ABSTRACT 

When  a  spherical  or  cylindrical  gas  cloud  expands  into  vacuum,  the 
resultant  flow  becomes  asymptotically  self-similar  and  independent  of 
the  original  cloud  diameter.  An  analytic  solution  for  the  asymptotic 
density  profiles,  given  by  Mirels  and  Mullen  [k]  is  re-examined  and  a  new 
solution  proposed,  which  gives  better  agreement  with  numerical  computations. 

Some  applications  of  the  asymptotic  solutions  are  discussed. 
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1 .  INTRODUCTION 

The  investigations  reported  here  were  initiated  by  a  paper  "Trans- 

# 

lational  Equilibration  of  Wires  Exploded  in  Vacua,"  by  0.  H.  Zihke  et  al  [l]  . 
The  authors  of  this  paper  made  use  of  an  unsteady  spherical  expansion  of 
ionized  gas  in  order  to  infer  the  temperature  of  the  gas  from  the 
characteristics  of  the  expansion.  The  method  described  by  Zinke  can  be 
a  powerful  tool,  provided  that  the  unsteady  flow  resulting  from  an  expanding 
gas  cloud  is  understood  theoretically.  It  is  believed  that  the  theoretical 
solution  presented  here  straightens  out  some  discrepancies  observed  by 
Zinke.  Further  remarks  on  this  are  given  in  section  5.  The  solution 
for  the  unsteady  cylindrical  expansion  can  also  be  applied  to  a  related 
problem,  the  expansion  of  a  hypersonic  jet  into  vacuum. 

In  the  above  mentioned  and  other  cases,  the  main  interest  is  in  the 
asymptotic  flow  field  after  long  times.  For  large  times  the  density 
distribution  of  the  expanding  gas  cloud  approaches  an  asymptotic  similarity 
form,  which  is  independent  of  the  original  diameter  of  the  cloud. 

It  can  be  shown,  however,  that  the  shape  of  the  asymptotic  density 
profile  depends  on  the  initial  density  distribution  of  the  gas  before 
expansion.  In  the  following  an  attempt  is  made  to  derive  approximate 
solutions  for  the  asymptotic  density  profiles  for  the  case  that  the  initial 
density  distribution  is  uniform.  Expansions  with  plane,  cylindrical  and 
spherical  symmetry  will  be  considered. 

1 - 

Numbers  in  brackets  denote  references  listed  on  page  25. 
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2.  ASYMPTOTIC  NATURE  OF  THE  FLOW 


i 


In  the  following  it  is  assumed  that  the  flow  is  isentropic  and  that 
the  gas  can  he  treated  as  a  continuum;  i.e.,  the  mean  free  path  is  very 
small  compared  with  the  initial  diameter  of  the  gas  cloud.  Then  the 
equations  of  motion  for  flows  of  cylindrical  (a  =  l)  or  spherical  (a  =  2) 


symmetry  are  respectively: 

§  =  o 

dt  dr  r 

(la) 

du  .  ■  du  .  1  dp 

aS  u  a?  is  a?  =  0 

(lb) 

(p/Po)  -  (o'/p,  )y 

(lc) 

p„  and  p0  are  taken  as  initial  density  and  pressure  if  the  density  is 
initially  uniform,  or  as  the  initial  values  at  r  =  0  if  the  density  is  not 
uniform.  It  is  assumed  that  the  gas  is  at  rest  initially  and  inclosed 
within  a  cylinder  or  sphere  of  radius  R0. 

For  the  following  we  introduce  nondimen sional  quantities: 


p  u  r  .  tan 

P  =  u=—  r  =  —  t  -  =-° 

Po  ®-0  “0 

where  ag  is  the  sound  speed  of  the  quiescent  gas  and  found  hy 

a  3  =  £o 
°  V  Po 


With  these  nondimensional  quantities  the  equations  of  motion  are: 

0 

•  0 


&+M  +  ,,E“ 

at  ar  r 


(2) 

(3) 

(Ua) 

(**> 


For  the  case  that  the  density  is  initially  uniform,  the  boundary  conditions 


are: 
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(5) 


at  t  =  0;  0  <  r  <  1:  p  =  1,  u  =  0 
1  <  r:  p  =  0 

at  t  £  0,  r  =  0:  u  =0 

For  the  case  that  a  =  0,  the  Equations  (4)  together  with  boundary  conditions 
(5)  describe  the  one -dimensional  (plane)  expansion  of  a  slab  of  gas 
extending  either  from  r  =  -1  to  r  =  +1  with  r  =  0  a  plane  of  symmetry  or 
from  r  =  0  to  r  =  1  with  a  fixed  wall  at  r  =  0.  In  a  few  cases,  exact 
solutions  of  Equations  (4a,  4b)  can  be  obtained.  These  solutions  which  we  shall 
discuss  in  the  next  paragraph,  will  serve  as  a  guide  to  construct 
approximate  asymptotic  solutions  for  the  more  general  case. 

Here  we  are  concerned  with  the  asymptotic  nature  of  the  flow  after  long 
times;  i.e,,  for  times  large  compared  to  R0/a0.  For  times  which  are  large 
in  this  sense,  the  original  gas  cloud  has  expanded  many  times  and  the 
pressure  term  pY~2  in  Equation  (4b)  becomes  relatively  small.  After 
neglecting  the  pressure  term  Equation  (4b)  is  solved  by 

u  =  r/t  (6) 

With  the  solution  (6)  for  u,  Equation  (4a)  has  the  solution: 

pta+1  *-  const  f(r/t)  (7) 

where  f(r/t)  is  ar  arbitrary  function  of  the  one  variable  r/t.  It  can  be 
seen  from  Equation  (6)  and  (7)  that  the  flow  becomes  asymptotically  self¬ 
similar;  i.e.,  the  density  and  velocity  distributions  at  different  times 
are  similar  and  differ  only  by  scaling  factors.  While  this  similarity  form 
of  Equations  (6)  and  (7)  means  a  certain  simplification  of  the  mathematical 
treatment,  the  problem  of  finding  the  particular  solution  for  f (r/t )  is  still 
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very  difficult  for  the  following  reason.  Since  any  solution  of  the  form 
(6)  and  (7)  solves  the  equations  of  motion  asymptotically,  any  of  these 
solutions  represents  a  possible  asymptotic  flow,  each  one  belonging  to 
different  initial  conditions.  Thus,  one  has  to  relate  a  given  initial 
density  distribution  to  a  corresponding  f(r/t)  and  it  seems  that  this 
can  be  done  exactly  only  by  computing  the  whole  flow  history  from  the 
beginning  of  the  expansion  up  to  large  times  where  the  flow  becomes  self- 
similar.  A  computation  of  this  sort  can  be  performed  numerically,  using  the 
finite  step  method  of  characteristics.  It  must  be  considered,  however, 
that  any  finite  step  computation  may  become  increasingly  inaccurate  when 
carried  out  for  too  many  steps,  due  to  the  accumulation  of  errors.  The 
number  of  steps  necessary  in  order  to  reach  the  asymptotic  similarity- 
solution  must  be  enormously  large.  In  spite  of  this,  numerical  computations 
for  uniform  initial  densities  and  various  values  of  y  and  o  have  been  carried 
out  by  several  authors  [2,3]  and  it  appears  that  the  results  of  these 
computations  are  surprisingly  accurate. 

3.  ANALYTICAL  SOLUTION 

A  first  attempt  toward  an  analytical  solution  has  been  made  by  Mirels 
and  Mullen  [4],  however  the  solution  does  not  show  satisfactory  agreement 
with  numerical  results  [2,3]. 

The  solution  proposed  here  is  obtained  by  only  slight  modification  of 
the  view  of  Mirels  and  Mullen  [4]  but  it  gives,  nevertheless,  far  better 
agreement  with  numerical  results.  The  procedure  used  to  derive  an  analytical 
expression  for  the  asymptotic  flow  consists  in  finding  a  valid  generalization 
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of  certain  exact  solutions.  As  has  "been  mentioned  before,  in  a  few  cases 
the  Equations  (4a),  (4b)  can  be  solved  exactly.  These  cases  are: 

1.  The  one -dimensional  (o  =  0)  expansion  of  an  initially  uniform 
slab  of  gas. 

2.  The  similarity  flows  for  c  =  0,  1,  2.  These  flows  are  strictly 
self-similar  at  all  times,  but  the  initial  density-distribution  is  not  uniform. 

The  characteristics  solution  for  one -dimensional  expansion  from  uniform 
initial  conditions  (case  l)  is  given  in  [5]  (Courant  and  Friedrichs)  arid 
is  discussed  in  greater  detail  in  [6,7).  In  the  asymptotic  limit,  this 
solution  simplifies  to: 


u  "  r/t  (8b) 

where  N  -  and  T  denotes  the  gamma-function. 

The  similarity  flows  have  been  investigated  in  References  [8]  and  [o]. 

If  R(t)  is  the  position  of  the  leading  edge  of  the  expansion,  one  can 

* 

define  a  similarity  variable 


y  =  r/R(t) 

so  that  the  gas  is  confined  in  0  s  y  <  1, 

Assuming  for  the  self-similar  flow  a  form  p  *  p"(a+1 ^f(y); 

v  -  H  .  g(y)  and  substituting  into  Equations  (4a,  b)  one  obtains 

l 

p  *  Rm^a+^(l  -  ya)Y“1' 

and  dR 

U  *  dt  *  y 


(9) 


(10) 

(11) 
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r 


i 


where  R  satisfies  the  differential  equation 


§  =  i  .  4  ti  -  r-^Xy-D]* 

dt  /d+1  y-1 

The  constants  of  integration  are  so  determined,  that  for 


t-°:  R=15t|=0 


In  the  asymptotic  limit  one  finds 


R  =  /TFT  *  *  * 

and  Equations  (10)  and  (ll)  become  asympt:-  Lically: 

pt”"1  .  (/Sn  t4)CT+1[i  -  (yorr  .  4  ff  f1 


U  =  V 


ptCT+1  =  d[~  1 


u  =  r/t 


Vvt  J  J 


(14a) 

(l4b) 


A  comparison  of  Equation  (8)  and  Equation  (l4)  shows  that  both  formulas 
are  quite  similar  and  of  the  form: 


(15a) 

(15b) 


where  D,  V  and  B  are  constants.  Thus,  we  assume,  in  accordance  with  Mirels 
and  Mullen  [4],  that  the  asymptotic  solution  for  uniform  initial  density  and 
ct  y  0  is  also  of  the  form  (15 ) ;  so  that  it  remains  only  to  determine  the 
constants  D,  V  and  B.  If  one  of  the  three  constants  is  given,  say  V, 
then  the  other  two  can  be  determined  from  the  conditions,  that  mass  and 
total  energy  are  conserved.  In  our  nondimen sional  variables  these  two  condi¬ 
tions  can  be  expressed  as: 
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„vt 

f.Pradr  =  ^-  ;  (16) 

0 

Jt 

.  o  I  “2  r° dr  =  Sw  •  •  (17) 

Equation  (l6)  equates  the  mass  of  the  expanded  gas  to  the  mass  at  t  =  0 
and  Equation  (17)  does  the  same  for  the  energy.  Note,  that  for  t  -»  «  the 
thermal  energy  vanishes  so  that  only  the  integral  over  the  kinetic  energy 
density  p/2  u2  is  taken  at  the  left  side  of  Equation  (17),  while  for  t  -  0 
the  kinetic  energy  is  zero  and  the  thermal  energy  per  unit  mass  is: 

Sq2 

CvT°  ~  Y(y-l) 

Introducing  the  new  variable  y  -  ^  and  inserting  (15a,  b)  into  (l6)  and  ( 17 ) 

gives: 

dvct+1  fa-y^Vay  -  ^  ;  (18) 

o 

DV^  f(l  -  jr*  )B  /+2dy  -  &  .  ^  .  (19) 

The  integrals  in  Equations  (l8)  and  (19)  can  be  transformed  into 
Eulerian  integrals  of  the  first  kind  and  expressed  in  terms  of  gamma 
functions.  It  is,  e.g., 

j  B  .  ros+un2^) 

f  ( l  "  y3 )  y  dy  -  4  ■■  . •  ( 20 ) 

Solving  Equation  (l8)  and  (19)  for  B  and  D  one  obtains: 
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!  2  r(B+1  +  IT) 

D  '  v"+1  '^'r(B+i)r(2i)  •  (22) 

Now  only  one  constant,  V,  remains  to  "be  determined  and  at  this  point  the 
present  analysis  differs  from  Mirels  and  Mullen  [4],  From  Equation  (15&) 
it  can  be  seen,  that  V  is  the  velocity  of  the  leading  edge  of  the 
expansion  (in  units  of  a0). 

Now  Mirels  and  Mullen  maintain  that  for  uniform  initial  density  the 

velocity  of  the  leading  edge  should  be  the  same  as  for  or.e-dimensional  flow., 

2 

and  therefore  V  =  >  regardless  of  the  value  of  a.  In  other  words,  the 

geometry,  whether  spherical  or  cylindrical,  does  not  influence  the  leading 

edge  velocity.  The  arguments  given  in  favor  of  this  are  difficult  and  not 

fully  convincing,  but  even  if  they  were  true  they  would  not  demand  for  V 

a  value  of  — j-  .  It  is  understood  that  Equation  (15a)  can  be  only  an 

approximate  formula  which  must  not  reflect  every  detail  of  an  exact  solution. 

It  could  be,  that  a  minute  fraction  of  the  total  mass  "escapes"  with  the 

2 

leading  edge  velocity  and  that  the  bulk  of  the  remaining  mass  expands 
with  an  "effective"  leading  edge  velocity  less  than  that.  In  order  to  have 
the  best  approximation,  valid  for  the  bulk  of  the  mass,  it  would  then  be 
advisable  to  pick  a  smaller  value  for  V.  This  is  also  suggested  by 
Equation  (l4a)  for  the  similarity  flow. 

Here  one  observes  that  the  geometry  does  have  an  influence  on  the 
leading  edge  velocity  expressed  by  the  geometrical  factor  •  If  the 
geometrical  factor  has  an  influence  in  the  case  of  similarity  flow, 
it  is  reasonable  to  assume  that  it  has  an  analogous  influence  in  the  case  of 
uniform  initial  density.  This  would  suggest: 
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V 


(23) 


1  2 

*  y-1  ’ 

for  both  similarity  and  uniform  initial  density  flow. 

A  more  rigorous  derivation  of  (23)  is  given  in  Appendix  A. 

With  (23)  inserted  into  (21)  and  (22)  the  constants  B  and  D  are: 


B 


)  -  ( 


D 


CT+1 


r(B+i  +  > 

r(B+l)  P(2|i)  ’ 


and  from  Equation  (15a), 

pt°+1  =  D  Cl  -  (£5  (Y-l) 


(2k) 

(25) 


(26) 


For  a  --  0  Equation  (8a)  is  recovered  from  Equation  (26);  i.e., 

Equation  (26)  is  exact  for  o  =  0.  For  a  j  0  the  formula  (26)  is  different 
from  Mirels  and  Mullen’s,  who  obtain  larger  values  for  the  exponent  B  and 
the  factor  D  as  a  result  of  the  larger  V  values  they  assumed.  In  fact, 
for  a  5  1  the  exponents B  obtained  by  Mirels  and  Mullen  are  even  larger  than 
the  exponent  for  similarity  flow,  which  is  (Equation  lUa).  A  consequence 
of  this  is  that  the  resultant  density  distribution  is  less  full  than  the 
corresponding  distribution  for  similarity  flow.  This  is  an  odd  result, 
considering  that  the  uniform  density  profile  is  fuller  initially  and  one 
ought  to  expect  that  it  also  remains  fuller  than  the  similarity  profile 
after  expansion.  The  exponents  obtained  by  Mirels  and  Mullen  are  so  large 
(tug.  for  o  *  2,  v  ==  7/5»  B  ■  8  instead  B  1  following  from  Equation  24) 
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that  they  essentially  cut  off  the  density  profile  at  a  distance  inside  of 
the  leading  edge  which  is  tantamount  to  our  supposition  that  only  a 
minute  fraction  of  the  mass  actually  achieves  a  velocity  close  to  that  of 
the  leading  edge.  In  other  words,  the  large  values  of  B  tend  to  partly 
compensate  for  the  prohibitively  large  leading  edge  velocity  assumed  by 
Mirels  and  Mullen. 

The  values  here  observed  for  B  (Equation  24)  are  always  smaller  and 
therefore  the  density  profiles  fuller,  than  for  the  corresponding  similarity 


flow. 

Apparently  the  formula  (26)  for  the  density  profile  becomes  nonsensical 
for  B  <  0,  or 


1  o+l 

Y-l  2 


(27) 


Wow,  it  is  interesting  to  note  that  inequality  (27)  describes  just  those  cases 

N+2 

which  are  physically  meaningless.  To  see  this,  we  recall  that  y  = 

where  N  is  the  number  of  internal  degrees  of  freedom  of  the  gas.  Therefore, 

On  the  other  hand  (o+l)  is  the  number  of  dimensions  of  the  space 
into  which  the  gas  expands.  If  this  is  (o+l)  =  n,  then  inequality  (27)  states 
N  <  n,  which  obviously  cannot  happen,  since  even  a  hypothetical  gas  must 
have  at  least  as  many  internal  as  external  degrees  of  freedom. 

The  density  factor  D  of  Equation  (25)  has  a  simple  physical  meaning: 
it  gives,  apart  from  the  time  factor  t  ,  the  density  at  the  origin  r  »  0. 

The  factor  D  is  plotted  in  Figure  1  and  Figure  2  for  o  =  1,  2  versus  y. 

For  comparison,  the  corresponding  curves  obtained  b  yMirels  and  Mullen  [4] 
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1.0  1.2  1.4  1.6  1.8  2.0 


FIGURE  1 

Density  factor  for  crsl 
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FIGURE  2 

Density  factor  for  <j  =  2. 
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and  the  results  of  numerical  computations  by  Greifinger  [2]  are  plotted 
in  the  same  diagram.  Also  the  curve  for  a  =  0  is  shown  in  Figure  1.  For 
this  case  the  results  of  Mirels  and  Mullen  [4],  Greifinger  [2]  and  the 
present  analysis  agree  with  each  other  and  the  exact  solution.  For 
a  =  1,  2  the  results  of  the  present  analysis  are  in  very  good  agreement 
with  the  numerical  data  [2]. 

Another  check  against  numerical  results  is  given  in  Figure  3,  which 
shows  the  density  profile  for  a  =  1,  y  =  5/3  in  comparison  with  the  results 
of  Mirels  and  Mullen  [4]  and  numerical  data  computed  at  NRL  [3].  All 
profiles  are  normalized  by  dividing  through  the  density  at  r  =  0.  Again, 
the  agreement  between  the  present  analysis  and  the  numerical  computation  is  very 
good.  In  Figure  2  the  values  for  D  according  to  Equation  (24)  are  plotted 
only  for  y  s  5/3.  For  larger  y-values  the  exponent  B  (Equation  24)  becomes 
negative  and  the  formula  (26)  is  no  longer  applicable.  One  might  suspect 
that  the  approximate  solution  (Equation  26)  becomes  less  accurate  on 
approaching  y  =  5/3. 

4.  LATERAL  EXPANSION  OF  A  HYPERSONIC  JET 
BOUNDED  BY  VACUUM 

If  the  lateral  flow  velocity  u  of  the  expanding  jet  is  small  compared 

with  its  axial  velocity  V0 ,  then  the  equations  of  motion  for  steady  flow 

reduce  to  Equations  (l)  with  t  replaced  by  ^  ,  where  x  is  the  axial 

vo 

distance  from  the  nozzle.  Using  the  nondimen sional  quantities  defined 
in  Equation  (2)  and  defining  a  nondimens ional  axial  distance  x  -  £  ,  one 
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numerical,  ref.(3) 


Asymptotic  density  profile  for 


Thus,  the  solution  for 


Xcl  X 

has  to  replace  t  in  Equations  (4)  hy  -rf-  =  r: 

"  o 

the  density  distribution  of  the  jet  is  obtained  from  Equation  (26)  for 
a  =  1  and  substituting  ^  for  t.  This  solution  is  valid  for  Mq  »  1. 

5.  MEASUREMENT  OF  PLASMA  TEMPERATURE 
In  their  paper.  "Translational  Equilibration  of  Wires  Exploded  in 

Vacua,"  [l]  the  authors  describe  a  method  for  measuring  the  temperature  of 
ionized  gas  with  the  aid  of  a  drift  tube.  The  method,  which  is  described 
in  detail  in  [1],  is  essentially  as  follows:  The  plasma  source,  in  this 
case  an  exploding  wire ,  is  located  at  one  end  of  a  long  evacuated  drift 
tube,  while  an  electrode  is  placed  at  its  opposite  end.  It  is  assumed 
that  the  wire  is  heated  uniformly  by  a  short  current  pulse.  After  the 
short  heating  period  the  gaseous,  ionized  wire  material  starts  to  expand 
without  being  heated  subsequently.  The  ion  flux,  after  reaching  the  electrode, 
gives  rise  to  a  time  dependent  current  in  the  electrode  circuit  which  can 
be  observed.  From  the  observed  current  distribution  the  temperature  of 
the  original  plasma  cloud  can  be  inferred. 

In  evaluating  the  experimental  data  Zinke  et  al.  [l]  assumed,  the 
plasma  would  expand  like  a  collision-free  gas;  i.e.,  according  to  free 
molecular  flow.  On  the  basis  of  this  assumption,  it  was  impossible  to  fit 
the  theoretical  prediction  to  the  experimental  data  unless  a  hypothetical 
"containment-time"  was  introduced.  In  other  word.:,  the  shape  of  the  observed 
current-time  distribution  was  quite  different  from  the  shape  of  any  of  the 
possible  theoretical  distributions.  The  theoretical  current  distributions 
for  free  molecular  flow  are  given  by 
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« 


(28) 


j(t‘)  =  fr«p(-^£) 

where  t  is  the  time  elapsed  from  the  explosion  of  the  wire,  r  the  distance 
between  the  explosion  and  the  electrode,  T  the  temperature  of  the  original 
plasma  cloud  and  K  a  constant  related  to  the  total  number  of  ions. 

The  distributions  Equation  (28)  form  a  two-parameter  class  of  curves 
with  parameters  K  and  T.  None  of  these  curves  fit  the  experimental  data. 

A  satisfactory  fit  could  be  achieved  only  after  introducing  a  third 
parameter,  H,  the  " containment - t ime Assuming  the  plasma  did  not  start 
to  expand  at  time  t  =  0,  i.e.,  immediately  after  heating,  but  at  a  later 
time,  t  =  H,  the  distribution  Equation  (28)  would  be  altered  to: 

d(t)  =  (!^hF  6XP[“'2M  fco®]  ‘  (29 

With  properly  chosen  K,  T,  H,  Zinke  could  obtain  a  reasonable  fit  to  all 

experimental  data.  However,  the  physical  meaning  of  the  containment-time, 

H,  remained  completely  unclear.  Obviously,  the  assumption  of  free  molecular 

flow  is  not  correct  for  the  high  initial  densities  encountered  in  the  wire 

explosion.  In  continuum  expansion  the  energy  of  the  gas  molecules  is 

redistributed  as  a  result  of  collisions.  This  may  significantly  change 

the  flow  and  the  current  distribution  observed  at  the  electrode  of  the 

drift  tube.  The  theoretical  current  distribution  for  continuum  flow  is 

proportional  to  pu  «  p  .  £  with  p  given  by  Equation  (26)  and  o  =  2.  Thus: 

t 


(30) 


Here  r  is  the  distance  between  the  electrode  and  the  plasma  source 

(exploding  wire)  and  t  the  time  elapsed  from  the  explosion. 

When  K  and  a0  were  chosen  properly,  a  good  fit  could  be  obtained  to  the 

experimental  data  of  Zinke  without  the  assumption  of  a  containment -time . 

kT 

Prom  the  temperature  can  be  determined  by:  =  y  —  . 

The  values  for  the  temperature  obtained  from  this  fit  were  about 
10s  °K.  These  values  are  more  likely  correct  than  those  obtained  by  Zinke 
of  about  106  CK. 


APPENDIX  A 

Introducing  a  new  dependent  variable  p*  by 

P  “  (P*)°+1  (31) 

and  substituting  (3l)  into  Equations  (4a,  b)  gives  the  following  set  of 
equations  for  p*: 

If  *  It  KM]  ■  0  <»> 

+  „  0  (32b) 

Let  us  now  consider  the  factor  (-  -  r^)  in  the  last  term  of  Eauation 

r  or 

(32a).  This  factor  is  exactly  zero  for  the  similarity  solution,  as  can  be 
seen  from  Equation  (ll). 

For  the  case  of  uniform  initial  density  this  facto*  becomes  asymptotically 
zero,  since  u  -  ■£.  It  is  also  zero  near  the  center  r  -  0,  which  follows 
from  the  boundary  conditions  (5)  according  to  which  u  •-  c.  r  +  cu  r2  +  . . . 

We  therefore  assume,  that  the  average  effect  of  the  last  term  of  Equation  (32a) 
is  small  and  to  a  certain  approximation  it  can  be  neglected  altogether. 
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With  the  further  transformations: 


=  jk  i  **  ■ 

(y*  -  2)  =  (y-2 ) (o+l )  +  a 

we  obtain  by  substituting  into  (32a ,b): 


(33) 


.  o  (34a) 

St*  Sr 

+  u*  +  (p*)<V*-2)  M.  ,  o  (34b) 

dt*  Sr  Sr 

From  (34a, b)  we  observe,  that  p*,  u*  satisfy  the  equations  of  motion  (4a, b) 
for  ct  =  0  with  y*  replacing  y.  p*,  u*  also  satisfy  the  same  boundary 
conditions  as  p,  u  for  the  case  of  uniform)  initial  density.  The  asymptotic 
solution,  therefore,  can  be  obtained  from  (8a, b)  with  p,  u,  t,  y  replaced 


by  p*,  u*,  t*,  y*5  i.e., 

p**  ~ 

(35a) 

u*  =  r/t* 

(35t) 

Taking  the  (o+l)  power  of  Equation  (35a)  and  substitution  of  (3l)  and 
(33)  into  (35a ,b)  gives: 


(3oa) 

(36b) 


The  constant  in  Equation  (36a)  can  now  be  re-adjusted  so  as  to  satisfy  the 
mass  conservation.  The  solution  (36a)  then  is  identical  with  Equation  (26). 
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